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Abstract 

We consider a DGP inspired brane scenario where the action on the brane is augmented 
by a function of the Ricci scalar, L{R). The cosmological implications that such a scenario 
\^ \ entails are examined for R n and shown to be consistent with a universe expanding with power-law 

acceleration. It is shown that two classes of solutions exist for the usual FRW metric and small 
Hubble radii. When the Hubble radius becomes larger, we either have a transition to a fully 5D 
regime or to a self-inflationary solution which produces a late accelerated expansion such that the 
radius becomes a function of n. 
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1 Introduction 
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The notion of extra dimensions and that they can be probed by gravitons and eventually non-standard 
matter has been the focus of attention in recent years. These models usually yield the correct 
Newtonian (l/r)-potential at large distances because the gravitational field is quenched on sub- 
millimeter transverse scales. This quenching appears either due to finite extension of the transverse 
dimensions [1, 2] or due to sub- millimeter transverse curvature scales induced by negative cosmological 
constants [3, 4, 5, 6, 7, 8]. A common feature of both of these types of models and also of the old 
Kaluza-Klein type models is the prediction of deviations from four-dimensional Einstein gravity at 
short distances. If the transverse length scale is not too small, this implies the possibility to generate 
bulk gravitons in accelerators or stars [2, 9]. The model of Dvali, Gabadadze and Porrati (DGP) 
[10], see [11] for extensions, is different, predicting that 4-dimensional Einstein gravity is a short- 
distance phenomenon with deviations showing up at large distances. The transition between four 
and higher-dimensional gravitational potentials in the DGP model arises as a consequence of the 
presence of both the brane and bulk Einstein terms in the action. Further, it was observed in [12, 13] 
that the DGP model allows for an embedding of the standard Friedmann cosmology in the sense 
that the cosmological evolution of the background metric on the brane can entirely be described by 
the standard Friedmann equation plus energy conservation on the brane. This was later generalized 
to arbitrary number of transverse dimensions in [14] . For a recent and comprehensive review of the 
phenomenology of DGP cosmology, the reader is referred to [15]. 

An interesting observation made a few years ago is that the expansion of our universe is currently 
undergoing a period of acceleration which is directly measured from the light-curves of several hundred 
type la supernovae [16, 17] and independently from observations of the cosmic microwave background 
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(CMB) by the WMAP satellite [18] and other CMB experiments [19]. However, the mechanism 
responsible for this acceleration is not well understood and many authors introduce a mysterious 
cosmic fluid, the so called dark energy, to explain this effect [20]. Recently, it has been shown that 
such an accelerated expansion could be the result of a modification to the Einstein-Hilbert action [21]. 
A scenario where the issue of cosmic acceleration in the framework of higher order theories of gravity 
in 4D is addressed can be found in [22]. One of the first proposals in this regards was suggested in 
[23] where a term of the form R~ l was added to the usual Einstein-Hilbert action. It was then shown 
that this term could give rise to accelerating solutions of the field equations without dark energy. 
In [24] a DGP brane model with a scalar field on the brane was proposed, predicting that for the 
gravitational potential, the mass density should be modified by the addition of the mass density of 
the scalar filed on the brane. Such a scalar field has its origin in the conformal transformation used 
to transform the action in the Jordan frame to the Einstein frame in the usual 4D modified gravity. 
For non-minimally coupled scalar field scenarios see [25]. 

In this paper, we focus attention on the DGP brane model where the action contains an arbi- 
trary function of the Ricci scalar, £(R), and obtain the evolution of the metric on the space-time. 
We concentrate on a specific form for C(R) = CqRJ 1 and solve the resulting dynamical equations, 
predicting a power-law acceleration on the brane. The components of the metric in the Gaussian 
normal coordinates are then calculated and presented. Finally, we show that there exists two classes 
of solutions close to the usual FRW cosmology for small Hubble radii in the model presented here. 
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2 DGP model with C(R) brane action 

We start by writing the action for the DGP model with an arbitrary function of the scalar curvature 
on the brane part of the action, that is 

J d 5 x^K + ^ J d 4 x^C(R) + S m [q^, Y> m ] , (1) 

where the first term in (1) corresponds to the Einstein-Hilbert action in 5D for the 5-dimensional 
bulk metric gAB, with the Ricci scalar denoted by 1Z. Similarly, the second term is the modified 
Einstein-Hilbert action with C{R) corresponding to the induced metric g M „ on the brane, where C{R) 
is some arbitrary function of the AD scalar curvature and 777,3 and 777,4 are reduced Planck masses in 
four and five dimensions respectively with S m being the matter action on the brane with a matter 
field denoted by ip m . The induced metric q^ is defined as usual from the bulk metric gAB by 

= tf$v9AB. (2) 

It would now be possible to write the field equations resulting from this action, yielding, in d — 1 
spatial dimensions 



7774 



TZab ~ \gABll) + ml« - ^uR^) S(y) = Sffi (f, u + T<£™>) 5(y), (3) 



C'{R) V T 

where T^ v = -pjj^T^ and T^ v is the energy-momentum tensor in the matter frame, and 

T^ urv) = {\q,u[C{R) - RC'(R)] + C'{R)^{q m q v p - q„ u q a p)} . (4) 

A prime here denotes differentiation with respect to R. Note that R lives on the brane. The cor- 
responding junction conditions, relating the extrinsic curvature to the energy-momentum tensor, 
become 



lim, 



_ C'{R)ml ( 1 a ^ R (d-i) 

m\ \ ^ 2(d-l) q ^ q U ^ 



(curv) 



y=0 

(5) 
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It should be noted that £{R) is a function of the brane parameters. 



3 Cosmology 

From the fact that the DGP model predicts deviations only at large distances, one might hope it 
could be ruled out from cosmological observations. However, as we shall see, it could account for 
cosmological equations of motion at any distance scale on the brane with any function of the Ricci 
scalar. Brane cosmology usually starts from the line element 



ds 2 = q^dx^dx" + b 2 {y,t)dy 2 = -n 2 (y,t)dt 2 + a 2 (y,t)-f ij dx i dx j + b 2 (y,t)dy 2 , 



(6) 



where -fij is a maximally symmetric 3-dimensional metric with k = —1,0,1 being the usual param- 
eters denoting the spatial curvatures. Building on the results of [26, 27], the cosmological evolution 
equations of a 3-brane in a 5D bulk resulting from equations (3) and (5) were presented in the first 
two references in [21]. Here we will follow [12] and give the results for a brane of dimension v + 1. 
Adopting the Gaussian normal system gauge 



b 2 (y,t) 



(7) 



The field equations on the brane for metric (6) and d = v + 1 spatial dimensions are 

,{v) 1 , .sofa 2 k 
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The junction conditions (5) for an ideal fluid on the brane, given by 

Tiu, = (p + p^Uv + pq^, 



n 2 a? 



(8) 
(9) 

(10) 



read 



where 
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pv-v = p + p^ 
p(tot) _ p + p (curv) 



(11) 

(12) 

(13) 
(14) 



and p and p are the energy density and pressure in the matter frame associated with Tu V respectively. 



Energy conservation on the brane follows from the vanishing of grfm Cos = =>■ v 
We obtain 



In' a af_ 

\n a a 



and in particular 



lim 

e^+0 



n 
n 

y=+e 
y=-e 



a' 



= lim 



r a 1 w 



y=-e 



= 0. 

(15) 
(16) 
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Insertion of (11) and (12) into this equation yields the equation of conservation 



p(tot) a 



y=0 



y=0 



(17) 



Also, insertion of (15) into the bulk equations Goo and G55 for y ^ yields a z/-dimensional version 
of the integral of [27], that is 



and 



1 2d 
-——Goo = =¥■ — a u Goo = t^I = 0, 
£ '(it) vn l ay 



— ^-G 55 = =t> -aa»G 55 = ~I = 0. 
£ (it) v at 



(18) 



(19) 



This means that if we define the quantities i + and i by taking the factor a v 1 out of the right hand 
sides of equations (18) and (19), that is 



/+= " -a' 2 + k a 



v-l 



-a' 2 + k\a 



u-l 



y>0 



y<0 



then I + and I are constants with the property that I + = I if 



lim a'\ v= + e = ± lim a'L=_ e . 



We can now simplify the previous equations by further restricting the gauge 

n(0,t) = l, 

and by simply performing the transformation 

t = J n{0,T)dr, 



(20) 
(21) 

(22) 
(23) 
(24) 



of the time coordinate. This gauge is convenient because it gives the usual cosmological time on the 
brane. Using equations (15) and (23), we find that our basic dynamical variable is a(y,t) with n(y,t) 
given by 

a(y,t) 



n(y,t) 



(25) 



d(0,t)" 

The basic set of cosmological equations in the present setting for any function of the Ricci scalar on 
the brane in the DGP model without a cosmological constant in the bulk are thus equations (12), (17), 
(20) and (21) which have to be amended with dispersion relations (or the corresponding evolution 
equations) for the ideal fluid components on the brane, that is 



+ 



k 



a(0,t) ' a(0,t) 



£{R){p + p( cu ™))a(0,t) 



vm 
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1+ = 
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a 2 (0,t)-a >2 (y,t) + k\ a^ l (y,t) 
a 2 (0,t)-a' 2 ( yi t) + k\ a v -\y,t) 



y>0 
y<0 ' 



(26) 

(27) 
(28) 
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(curv) 



C'{R) 



-[£(R)-R£'(R)}-3 



d(0,t) 
a(0,t) 



RC"(R) , 



(29) 



P (curv) = i^{ 2 (^|) *£"( R ) + + R 2 £"'( R ) - - > (30) 

a(y,t) 



n(y,t) 



(31) 



d(0,t)" 

Let us now discuss the cosmology in the DGP model by taking I + = I~ . The cosmological equations 
in this framework for a (v — l)-dimensional space are given by 



a 2 (0,t) + k 2(p + p( cu ™)) 



a 2 {0,t) 



v(y — l)rriv 



v-l ' 



a 2 (0,t)-a r \y,t) + k 
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(y,t), 



n(y,t) 



d(0,t)' 



(32) 

(33) 
(34) 



The evolution of the background geometry of the observable universe according to the Friedmann 
equation can thus be embedded in the DGP model, with the behavior of a(y, t) off the brane de- 
termined solely by the integral I and the boundary condition a(0, t) from the Friedmann equation. 
This embedding will be asymmetric in all realistic cases, because the requirement that the Friedmann 
equation holds on the brane is equivalent to the smoothness condition 



lim a'(e,i) = lim a'(—e,t). 



(35) 



This could yield a symmetric embedding only for a'(0,t) = 0, but this is incompatible with the time 
independence of the integral I apart from the case k = — 1, a = 1. For a discussion of this point the 
reader may consult [28, 29]. In the present calculations we will choose the sign of y in the direction 
of the increasing scale factor, a' > 0. The possibility of a direct embedding of Friedmann cosmology 
is a consequence of the fact that the evolution of the background geometry (6) and the source terms 
p(curv) anc j p(curv) are SU pp 0sec j to depend only on t and y. For v = 3 we obtain 



d 2 (0,t) -a f2 (y,t) + k a z (y,t), 



(36) 



and from the equation for n(y,t), the solutions for the metric components off the brane in terms of 
the metric on the brane (assuming a' > 0) are 



a 2 (y, t) = a 2 (0, t) + (d 2 (0, t) + k)y l + 2 v /(d 2 (0, t) + k)a 2 (0, t) - Iy, 



n(y,t) 



a(0, t) +d(0, t)y 2 + a(0,t)y 



a(Q,t)a(0,t) +d 2 (0,t) + k 



(37) 



(38) 



^(a^(0,t) + k)a 2 (0,t)-l\ a(y,ty 

This embedding of the Friedmann cosmology on the brane becomes particularly simple for 1 = 0, 
that is 

(39) 



a{y,t) = a(0,t) + yja 2 (0,t) + ky, 
d{0,t) 



n(y,t) = 1 + 



v/d 2 (0,t) + k 



y- 



(40) 



5 



4 R n gravity in the DGP model 

To progress further, the form of C(R) should be specified. For ease of exposition and clarity, let us 
focus attention on theories where a R n term is present in the action and write 

C{R) = C Q R n . (41) 

Let us also assume a power law solution form for the brane scale factor as 

a(0,t) = a ^y. (42) 

The interesting cases are for a > 1 which give rise to acceleration. 

To proceed, we consider the evolution of the scale factor with time on the brane. Using equation 
(32) for the spatially flat FRW metric and setting v = 3, we write 



2 

Now, use of the conservation equation leads to 



- 1 ,<**>. (43) 



a(0,t) 1 



a(0, t) 6m| 



( tot ) _|_ 3p(* oi ) . (44) 



We must now solve the system of equations (43) and (44) with p = p = 0. Substituting equations 
(41) and (42) into the above equations describing the dynamical system for the brane, we obtain an 
algebraic system for the parameters n and a 

a[a(n - 2) + 2n 2 - 3n + 1] = 0, 
a[a(n - 2) + (n 2 - n + 1)] = n(2n - l)(n - 1), (45) 



from which the allowed solutions 



a = => n = 0, 1/2, 1, 



2n 2 - 3n + 1 

a = , V n except n = 2, (46) 

2 — n 

follow. The solutions with a = are not interesting since they provide static cosmologies with a 
non-evolving scale factor on the brane, matching cosmological models resulting from the solutions 
of the Einstein equations without matter and n = 1. This particular value of n warrants a further 
discussion. It is well known that in ordinary 4D gravity an action of the form C{R) = CqR 71 results in 
a constant scale factor in a cosmological setting without any ordinary matter and a singular equation 
of state for n = 1 [22]. It is therefore not surprising to expect the same behavior in DGP models when 
the action on the brane is taken as that mentioned above. However, if one takes C(R) = R + CoR n , 
one is lead to ordinary general relativity as the low energy limit of the theory and one finds that 
this theory is equivalent to scalar quintessence models through conformal transformations [23]. The 
corresponding result for C(R) = CoR n , n > and n ^ 1 is also spelled out in [22]. Using equations 
(29) and (30) we can deduce the equation of state for the family of solutions a / 0. We then have 

{curv) _ ( 6n 2 -7n-l\ 

which clearly behaves as w^ curv ^ — > — 1 for n — > oo, playing the role of a cosmological constant. For 
R~ l {n = —1) which corresponds to a = 2 we have a powerdaw acceleration on the brane without 
having to introduce dark energy. This result is consistent with the observational results similar to 
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Figure 1: Behavior of q(n) as a function of n. 

dark energy with the equation of state parameter —1 < w^ curv ^ = < — | [30]. Using equations 
(42) and (46) we obtain the deceleration parameter 

, /2n 2 -2n-l\ 
g(n) = -( 2 n2-3n + l J- (48) 

It is clear that for n — ► oo we have g(oo) = —1. This means that the universe (brane) is continuously 
expanding at an ever- increasing rate. Clearly, for n = — 1 or a = 2 we obtain q = — \. It is worth 
noting again, as was mentioned before, that in the normal DGP model, that is when n = 1 and 
with no ordinary matter present, we cannot define the equation of state and deceleration parameter 
because the scale factor on the brane is constant, see equations (43) and (44). We therefore expect 
the same behavior in our model and see that these equations diverge for n = 1. Figure 1 shows the 
behavior of q as a function of n. As can be seen, for n < —0.4 and n > 1.4, q — > — 1. 

Let us now find the evolution of a(y, t) and n(y, t) everywhere in space-time. Substituting equa- 
tions (42) and (46) into equations (37) and (38), one finds 




The scale factor on the full space-time, a(y, t), is plotted in figure 2 for the special case 1 = and a = 2 
or n = —1. Note that for y = 0, equations (51) and (52) reduce to (42) and n(0,t) = 1 respectively. 
There appears coordinate singularities on the space- like hypercone y = ±(^j). This is presumably 
a consequence of the fact that the orthogonal geodesies emerging from the brane (which we used to 
set up our Gaussian normal system, b 2 = 1) do not cover the full five-dimensional space-time. 
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Figure 2: Behavior of a(y,t) as a function of t and y for n = — 1. 



5 Brane Friedman equations with C(R) 

In this section we derive the Friedman equations for the brane metric in the presence of C(R) in the 
brane part of DGP action. Thus, for I = in equation (36) and using equation (23) we can write 

a 2 (0,t)-a' 2 (y,t)}a 2 (y,t) = 0. (53) 

Note that we assume k = and the bulk cosmological constant is zero. Taking the symmetry y <-> — y 
(I + = I~) for simplicity, equation (26) can be used to compute a'(y, t) on the two sides of the brane. 
We have in this case [a'(y, t)] = 2a'(0 + ,t). By continuity when y — > 0, equation (53) yields the 
generalized (first) Friedman equation 

#2_ 2e _V~ = !*L p {tot)^ (54) 



k 2 C'(R) 3 
M K 2 _ m -: 

a(0,t)n(0,t) ' K ~ " L 5 ' 



where the Hubble parameter is H = 7Kfr% , k 2 = m 5 3 , = m 4 2 and e = ±1 is the sign of 



«'(//•')• 

Now, let us discuss the solutions of the Friedman equation (54) together with (17) when the bulk 
cosmological constant vanishes. It is now apparent from (54) that the standard cosmology, namely 
the usual AD Friedman equation (43) is recovered whenever the last term on the left hand side of 
(54) is subdominant with respect to the first term, namely when 

H>>2 i?m- (55) 

Thus using ansatzs (41) and (42), the above equation in terms of the Hubble radius H -1 can be 
written 

, 1 n-l 

H- 1 < (r c n£ ) — (12 - 6/a)~, (56) 

where r c = t^V- For n = 1 this matches the scale r c found in [10], setting the crossover between the 
AD and 5D gravity regimes. For the special case n = — 1, ora = 2 for which a power-law acceleration 
consistent with t 2 is obtained we have 

H- 1 « (81)3(-r c £ )^. (57) 
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Also, for n — ► 2 (a — ► oo) one has 

i^- 1 < (24r c £ )^ (58) 

Thus, the Hubble radius if -1 in this model depends on n. For the Hubble radius when equation 
(55) is not satisfied, there are two distinct behaviors depending on e. Equation (54) can indeed be 
rewritten as 

"-^^i^'^' (59) 

Let us start by examining the case where e = —1. Using condition HC'(R) <C r" 1 , one may then 
expand (54) to obtain 



H K 2 



C'{R) 6 



-p^ tot \ (60) 



or 

&C'{R)-* = ^pWr, (61) 



which is the full 5D regime with C{R) gravity on the brane. Equation (61) for R n becomes 

^ 4 " +6 = ^(n£ ) 2 (l2- (p + p^f. (62) 



Note that p( curv ) is a function of R. One has thus a transition from a 4D to a 5.D regime. If e = 1, 
however, H is always larger than -ff se if given by 



2u 2 

flself = "4, (63) 



where # se if = HC'(R). For i?" we find 



/ 2/i 2 \ l-n 

and the expansion will never enter into a fully 5D regime. Thus for fc = or fc = —1 the Hubbel 
parameter is bounded from below by equation (64) for R n on the brane. 

An interesting feature of the normal DGP models is the existence of ghost-like excitations [31, 
32, 33]. It would therefore be interesting to briefly mention recent results relevant to the present 
work. In [33], the author has shown, in the context of normal DGP models, that if we introduce a 
positive cosmological constant on the brane (Hr c > 1), then the spin-2 graviton will have a mass 
in the range < m 2 < 2H 2 and that there is a normalizable brane fluctuation mode with mass 
m 2 = 2H 2 . Although the brane fluctuation mode is healthy, the spin-2 graviton has a helicity-0 
excitation which is known as a ghost. If we allow a negative cosmological constant on the brane, the 
brane bending mode becomes a ghost for 1/2 < Hr c < 1. This confirms the results obtained by the 
boundary effective action that there exists a scalar ghost mode for Hr c > 1/2. In a self-accelerating 
universeiir c = 1, the spin-2 graviton has mass m 2 = 2H 2 , which is known to be a special case for 
massive gravitons in de Sitter spacetime where the graviton has no helicity-0 excitation and so no 
ghost. However, in DGP models, there exists a brane fluctuation mode with the same mass and 
there arises a mixing between the brane fluctuation mode and the spin-2 graviton. Thus, according 
to this scenario, our model may be prone to having ghosts and this could be the subject of a separate 
investigation. 



9 



6 Conclusions 



Brane models generally provide an interesting extension of our parameter space for gravitational 
theories. In this work we have discussed the DGP model with an arbitrary function of the Ricci 
scalar in the brane part of action. The cosmological evolution of this model was studied by solving 
the relevant dynamical equations. The components of the metric on the space-time was obtained 
for the R n term in 4D gravity. The evolution of the universe in such a scenario was shown to be 
consistent with the present observations, predicting an accelerated expansion. Finally we have shown 
that there exists two classes of solutions, close to the usual FRW cosmology for small enough Hubble 
radii for C(R) gravity. 



References 

[I] N. Arkani-Hamed, S. Dimopoulos, G. Dvali, Phys. Lett. B 429 (1998) 263, [hep-ph/9803315], 
N. Arkani-Hamed, S. Dimopoulos, G. Dvali, Phys. Rev. D 59 (1999) 086004, [hep-th/9807344]. 

[2] I. Antoniadis, N. Arkani-Hamed, S. Dimopoulos, G. Dvali, Phys. Lett. B 436 (1998) 257, [hep- 
ph/9804398]. 

[3] L. Randall, R. Sundrum, Phys. Rev. Lett. 83 (1999) 4690, [hep-th/0906064]. 

[4] W. Muck, K. S. Viswanathan, I. V. Volovich, Phys. Rev. D 62 (2000) 105019, [hep-th/0004017]. 

[5] R. Gregory, V. A. Rubakov, S. M. Sibiryakov, Class. Quan. Grav. 17 (2000) 4437, [hep- 
th/0003109]. 

[6] I. Ya. Aref'eva, M. G. Ivanov, W. Muck, K. S. Viswanathan, I. V. Volovich, Nucl. Phys. B 590 
(2000) 273, [hep-th/0004114]. 

[7] M. Cvetic, M. J. Duff, J.T. Liu, H. Lu, C. N. Pope, K. S. Stelle, Nucl. Phys. B 605 (2001) 141, 
[peh-th/0011167]. 

[8] B. Abdesselam, N. Mohammedi, Phys. Rev. D 65 (2002) 084018, [hep-th/0110143]. 

[9] E. A. Mirabelli, M. Perelstein, M. E. Peskin, Phys. Rev. Lett. 82 (1999) 2236, [hep-ph/981133], 
T. Appelquist, H. C. Cheng, B. A. Dobrescu, Phys. Rev. D 64 (2001) 035002, [hep-ph/0201131], 
T. G. Rizzo, Phys. Rev. D 64 (2001) 095010, [hep-ph/0106336], 
S. Cullen, M. Perelstein, Phys. Rev. Lett. 83 (1999) 268, [hep-ph/9903422], 
V. Barger, T. Han, C. Kao, R. J. Zhang, Phys. Lett. B 461 (1999) 34, [hep-ph/9905474], 
S. Cassisi, V. Castellani, S. DegPInnocenti, G. Fiorentini, B. Ricci, Phys. Lett. B 481 (2000) 323, 
[astro-ph/0002182], 

M. Biesiada, B. Malec, Phys. Rev. D 65 (2002) 043008, [astro-ph/0109545], 

S. Hannestad, G. G. Raffelt, Phys. Rev. Lett. 88 (2002) 071301, [hep-ph/0110067]. 

[10] G. Dvali, G. Gabadadze, M. Porrati, Phys. Lett. B 485 (2000) 208, [hep-th/0005016]. 

[II] G. Dvali, G. Gabadadze, Phys. Rev. D 63 (2001) 065007, [hep-th/0008054], 

G. Dvali, G. Gabadadze, M. Kolanovic, F. Nitti, Phys. Rev. D 65 (2002) 024031, [hep-th/0106058]. 

[12] R. Dick, Class. Quant. Grav. 18 (2001) Rl, [hep-th/0105320]. 

[13] R. Dick, Actaphys. Polon. B 32 (2001) 3669, [hep-th/0110162]. 

[14] R. Cordero, A. Vilenkin, Phys. Rev. D 65 (2002) 083519, [hep-th/0107175]. 

[15] A. Lue, Phys. Rept. 423 (2006) 1, [astro-ph/0510068]. 



10 



[16] A. G. Riess et. al [Supernova Search Team Collaboration], Astron. J. 116 (1998) 1006, [astro- 
ph/9805201]. 

[17] S. Perlmutter et. al, Astron. J. 517 (1999) 565, [astro-ph/9812133], 

D. N. Spergel et. al, Astrophys. J. Suppl. 148 (2003) 175, [astro-ph/0302209]. 

[18] C. L. Bennett et. al, Astrophys. J. Suppl. 148 (2003) 1, [astro-ph/0302207]. 

[19] C. B. Netterfiled et. al, Astrophys. J. 571 (2002) 604, [astro-ph/0104460], 
N. W. Halverson et. al, Astrophys. J. 568 (2002) 38, [astro-ph/0104489]. 

[20] S. M. Carroll, Living Rev. Rel. 4 (2001) 1, [astro-ph/0004075]. 

[21] C. Deffayet, Phys. Lett. B 502 (2001) 199, [hep-th/0010186], 

C. Deffayet, G. R. Dvali and G. Gabadadze, Phys. Rev. D 65 (2002) 044023, [astro-ph/0105068], 

C. Deffayet and S. J. Landau, J. Raux, M. Zaldarriaga and P. Astier, Phys. Rev. D 66 (2002) 
024019, [astro-ph/0201164], 

J. S. Alcaniz, Phys. Rev. D 65 (2002) 123514, [astro-ph/0202492], 

D. Jain, A. Dev and J. S. Alcaniz, Phys. Rev. D 66 (2002) 083511, [astro-ph/0206224], 

A. Lue, R. Scoccimarro, G. Starkman, Phys. Rev. D 69 (2004) 044005, [astro-ph/0307034]. 

[22] S. Capozziello, Int. J. Mod. Phys. D 11 (2002) 483, [astro-ph/0201033], 
S. Capozziello, Int. J. Mod. Phys. D 12 (2003) 1969, [astro-ph/0307018]. 

[23] S. M. Carroll, V. Duvvuri, M. Trodden, M. Turner, Phys. Rev. D 70 (2004) 043528, [astro- 
ph/0306438]. 

[24] K. Atazadeh and H. R. Sepangi, Phys. Lett. B 643 (2006) 76, [gr-qc/0610107]. 

[25] K. Nozari, Phys. Lett. B 652 (2007) 159, [hep-th/0707.0719], 
K. Nozari, JCAP 09 (2007) 003, [hep-th/0708.1611]. 

[26] P. Binetruy, C. Deffayet, D. Langlois, Nucl. Phys. B 565 (2000) 269, [hep-th/9905012]. 

[27] P. Binetruy, C. Deffayet, U. Ellwanger, D. Langlois, Phys. Lett. B 477 (2000) 285, [hep- 
th/9910219]. 

[28] P. P. Avelino, C. J. A. P. Martins, Astrophys. J. 565 (2002) 661, [astro-ph/0106274]. 

[29] C. Deffayet, G. Dvali, G. Gabadadze, [astro-ph/0 106449]. 

[30] T. Padmanabhan, Phys. Rept. 380 (2003) 235, [hep-th/ 0212290]. 

[31] M. A. Luty, M. Porrati and R. Rattazzi, JHEP 0309 (2003) 029, [hep-th/ 0303116]. 

[32] A. Nicolis and R. Rattazzi, JHEP 0406 (2004) 059, [hep-th/0404159]. 

[33] K. Koyama, Phys.Rev. D 72 (2005) 123511, [hep-th/0503191]. 



11 



